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ABSTRACT 

We assess supernova (SN)-driven pregalactic outflows as a mechanism for distributing the product 
of stellar nucleosynthesis over large cosmological volumes prior to the reionization epoch. SN ejecta 
will escape the grasp of halos with virial temperatures T vlr > 10 43 K (corresponding to masses M > 
10 8 hr 1 M© at redshift z — 9 when they collapse from 2-a fluctuations) if rapid cooling can take place, 
and a significant fraction of their baryonic mass is converted into stars over a dynamical timescale. We 
study the evolution of SN-driven bubbles as they blow out from subgalactic halos and propagate into 
the intergalactic medium (IGM) , and show that to lift the halo gas out of the potential well the energy 
injection must continue at least until blow-away occurs. If the fraction of ionizing photons that escape 
the dense sites of star formation into the intergalactic space is greater than a few percent, pregalactic 
outflows will propagate into an IGM which has been pre-photoionized by the same massive stars which 
later explode as SNe, and the expansion of the metal-enriched bubbles will be halted by the combined 
action of external pressure, gravity, and radiative losses. The collective explosive output of about ten 
thousands SNe per M > 10 8 h~ x M© halo at these early epochs could pollute vast regions of intergalactic 
space to a mean metallicity (Z) = fiz/flj, > 0.003 (comparable to the levels observed in the Lya forest 
at z « 3) without much perturbing the IGM hydrodynamically, i.e. producing large variations of the 
baryons relative to the dark matter. Rayleigh- Taylor instabilities between the dense shell which contains 
pristine swept-up material and the hot, metal-enriched low-density bubble may contribute to the mixing 
and diffusion of heavy elements. The volume filling factor of the ejecta is higher than 20% if the star 
formation efficiency is of order 10%. Larger filling factors (not required by current observations) may 
be obtained for larger efficiencies, moderately top-heavy IMFs, halos where a significant fraction of the 
gas is in a galactic disk and does not couple to the outflow (as matter is ejected perpendicularly to the 
disk), or from a population of more numerous sources - which would therefore have to originate from 
lower amplitude peaks. When the filling factor of the ejecta becomes significant, enriched material will 
typically be at a higher adiabat than expected from photoionization. 

Subject headings: cosmology: theory - galaxies: formation - intergalactic medium quasars: 
absorption lines 



1. INTRODUCTION 

In currently popular cosmological scenarios - all variants 
of the cold dark matter (CDM) cosmogony with differ- 
ent choices for the parameters Qmi ^a, h, as, n - some 
time beyond a redshift of 15 the gas within halos with 
virial temperatures T vu > 10 4 K [or, equivalently, M > 
10 9 (1 + z)~ 3 / 2 /i~ 1 M Q ] 5 cooled rapidly due to the excita- 
tion of hydrogen Lya by the Maxwellian tail of the electron 
distribution, and fragmented. Massive stars formed with 
some initial mass function (IMF), synthesized heavy ele- 
ments, and exploded as Type II supernovae (SNe) after a 
fewxf0 7 yr, enriching the surrounding medium; these sub- 
galactic stellar systems, aided perhaps by an early popula- 
tion of accreting black holes in their nuclei, generated the 
ultraviolet radiation and mechanical energy that reheated 
and reionized the universe. While collisional excitation 
of molecular hydrogen may have allowed the gas in even 
smaller systems [virial temperatures of only a few hundred 



K, corresponding to masses around 10 7 (I+z) 3/>2 h 1 M©] 
to cool and form stars at earlier times (Couchman & Rees 
f986; Haiman, Rees, & Loeb 1996; Tegmark et al. 1997), 
H 2 molecules are efficiently photo-dissociated by stellar 
UV radiation, and such 'negative feedback' could have 
suppressed molecular cooling and further star formation 
inside very small halos (e.g. Haiman, Abel, & Rees 2000; 
Ciardi, Ferrara, & Abel 2000). 

Throughout these crucial formative stages, the all- 
pervading intergalactic medium (IGM) acted as a source 
for the gas that is accreted, cools, and forms stars within 
these subgalactic systems, and as a sink for the metal- 
enriched material, energy, and radiation which they eject. 
The well-established existence of heavy elements like car- 
bon, nitrogen, and silicon in the Lya forest clouds at 
z — 3 — 3.5 may be the best evidence for such an early 
episode of pregalactic star formation. The detection of 
weak but measurable C iv and Si iv absorption lines in 
clouds with H i column densities as low as I0 14 5 cm -2 
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implies a minimum universal metallicity relative to solar 
in the range [—3.2] to [—2.5] at z = 3 — 3.5 (Songaila 
1997). There is no indication in the data of a turnover 
in the C iv column density distribution down to A^iv ~ 
10 117 cm" 2 (Afci « 10 142 cm" 2 , Ellison et al. 2000); the 
analysis of individual pixel optical depths may actually im- 
ply the presence of weak metal lines below the detection 
threshold (Cowie & Songaila 1998; Ellison et al. 2000). 
Widespread enrichment is further supported by the recent 
observations of O vi absorption in low density regions of 
the IGM (Schaye et al. 2000). 

In this paper we argue that the observed paucity of re- 
gions in the IGM which are of truly primordial compo- 
sition - or have abundances as low as those of the most 
metal-poor stars in the Wilky Way halo (Ryan, Norris, 
& Beers 1996) - may point to an early enrichment epoch 
by low-mass subgalactic systems, rather than being due 
to late pollution by massive galaxies. Whilst outflows 
of metal-rich gas are direcly observed in local starbursts 
(e.g. Heckman 1999) and z w 3 Lyman-break galaxies 
(Pettini et al. 2000) , most of this gas may not leave these 
massive galaxies altogether, but remains trapped in their 
gravitational potential wells until it cools and rains back 
onto the galaxies. By contrast, metal-enriched material 
from SN ejecta is far more easily accelerated to veloci- 
ties larger than the escape speed - about 80 kms -1 at 
the center of a lO 8 ^ 1 M Navarro-Frenk- White (1997, 
hereafter NFW) halo at z — 9 - in the shallow potential 
wells of subgalactic systems (Larson 1974; Dekel & Silk 
1986). These early protogalaxies - with masses compa- 
rable to those of present-day dwarf ellipticals - are then 
expected to release significant amounts of kinetic energy, 
heat, and heavy elements into the surrounding intergalac- 
tic gas. Many authors have addressed the impact of SN- 
driven winds from small starbursting galaxies at high-z 
(e.g. Barkana & Loeb 2001; Ferrara, Pettini, & Shchekinov 
2000; Scannapieco & Broadhurst 2001; MacLow & Ferrara 
1999; Murakami & Babul 1999; Nath & Trentham 1997; 
Voit 1996; Tegmark, Silk, & Evrard 1993) or from massive 
galaxies at more modest redshifts (Aguirre et al. 2000; 
Theuns, Mo, & Schaye 2001) on the thermal and chemical 
state of the IGM. Gnedin & Ostriker (1997) and Gnedin 

(1998) have argued that violent merging between proto- 
galaxies may be an alternative mechanism to transport 
metals in the low density IGM, while Barkana & Loeb 

(1999) have discussed the possibility that enriched ma- 
terial in halos with virial temperatures < 10 4 K could be 
photoevaporated at reionization. Here we give an idealized 
assessement of pregalactic outflows as a mechanism for dis- 
tributing the products of stellar nucleosynthesis over large 
volumes. We will be concerned exclusively with explosive 
multi-SN events operating on the characteristic timescale 
of a few x 10 7 yr, the lifetime of massive stars. In a com- 
plementary study, Efstathiou (2000) has recently shown 
that a large fraction of the baryonic mass of galaxies with 
virial temperatures ss 10 5 K can also be expelled in 'quies- 
cent' mode, i.e. over the relatively long timescale of 1 Gyr. 
To anticipate the conclusions of this work, we find that the 
IGM will be polluted over large scales at the end of the 
'dark ages' if a fraction greater than a few percent of the 
baryonic mass of subgalactic halos can be converted into 

6 Throughout this work we adopt the notation Y x = Y/10 x , and eg! 



stars over a dynamical timescale. In the case of large star 
formation efficiencies and/or moderately top-heavy IMFs, 
the temperature of vast regions of the IGM will be driven 
to a higher adiabat than expected from photoionization, 
so as to inhibit in these regions the formation of further 
protogalaxies by raising the Jeans mass. 

2. AN EARLY ENRICHMENT EPOCH? 

Before embarking on a discussion of the consequences of 
pregalactic outflows on the thermal and chemical state of 
the IGM, it is worth summarizing a few key observational 
facts and theoretical results: 

(1) Numerical N-body/hydrodynamics simulations of 
structure formation in the IGM within the framework of 
CDM dominated cosmologies (e.g. Cen et al. 1994; Zhang, 
Anninos, & Norman 1995; Hernquist et al. 1996; Theuns 
et al. 1998) have recently provided a coherent picture for 
the origin of the Lya forest, one of an interconnected net- 
work of sheets and filaments with virialized systems (ha- 
los) located at their points of intersection. The simula- 
tions show good agreement with the observed line statistics 
under the assumption that an IGM with baryon density 
parameter Qbti 2 = 0.019 (Buries and Tytler 1998) is pho- 
toionized and photoheated by a UV background with a hy- 
drogen ionization rate T(z = 3.5) « 0.6 x 10~ 12 s _1 , close 
to that inferred from quasars (Haardt & Madau 1996). 
They also show a clear correlation between H i column, 
gas temperature, and overdensity, with 

p b /p b ~ 0.8 JV££ 13 and p b /p b ~ 0.3 T 2 (1) 

at z = 3.5 (Ricotti, Gnedin, & Shull 2000; Zhang et al. 
1998; Hui & Gnedin 1997). 6 In photoionization equilib- 
rium, an optically thin cloud with internal density p b will 
have a neutral hydrogen fraction of 

where the temperature dependence is from the radia- 
tive recombination rate. Combining equations (|l]) and 
(|^), and omitting for simplicity the redshift and pho- 
toionization rate scalings, one obtains approximately 

^hi/^H — 10 ~ 5 ' 3 A^j 2 13 . The baryonic mass fraction in 
the Lya forest per unit logarithmic H i column interval 
can be written as 

_ 1.3m p N 2 J{N m ,z) dz_ 
"hPcrit nm/nn cat 

where m p is the proton mass, p cr jt = 3H 2 /8irG is the 
critical density of the universe at redshift z, f(Nm,z) — 
10~ 121 (1 + z ) 2 ' 5 is the bivariate distribution of 

H i columns and redshifts (cf. Kim et al. 1997), and cdt/dz 
is the line element in a Friedmann cosmology. From these 
relations one finds a mass fraction f m ~ 0.6 h which is 
independent of H i column, i.e. most of the baryons at 
this epoch lie in the range 12.5 < logA^i < 14.5 and are 
distributed equally per decade in column density (Zhang 
et al. 98). The volume filling factor of Lya forest clouds is 

fv = f m ^~ 0.75 hN^l (4) 

Pb 

units. 



Madau, Ferrara & Rees 



3 



Hence, in this picture, the metals associated with 
logN ffl < 14.2 filaments fill a fraction > 3% (h = 0.5) 
of intergalactic space, and are therefore far away from 
the high overdensity peaks where galaxies form, gas cools, 
and star formation takes place. Their chemical enrichment 
must then reflect more uniform (i.e. 'early') rather than 
in-situ (i.e. 'late') metal pollution. 

(2) At z = 3 - 3.5, clouds with N m Z 10 147 cm" 2 show a 
spread of at most an order of magnitude in their metallic- 
ity, and their narrow line widths require that they be pho- 
toionizcd and cold rather than collisionally ionized and hot 
(Songaila & Cowie 1996). At these redshifts, hot rarefied 
gas, exposed to a metagalactic ionizing flux, will not be 
able to radiatively cool within a Hubble time. The cool- 
ing time can be defined as the ratio of the specific energy 
content to the radiative cooling rate, 



could be explained if the ejection of heavy elements at 
velocities exceeding the small escape speed of subgalactic 
systems were to take place at very high redshifts. Hot en- 
riched material cools more efficiently at these early epochs, 
since Ht coo \ oc (1 + z)~ 3 / 2 and the Compton cooling time 
of the shocked ionized ejecta off cosmic microwave back- 
ground (CMB) photons, 
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where k is Boltzmann's constant, T is the gas tempera- 
ture, and A is the radiative cooling function. Here we have 
computed A for a primordial plasma 7 with helium frac- 
tion by mass equal to 0.25, and total number density of all 
species n. The cooling function, based on the rates given 
in Hui & Gnedin (1997), was calculated in the presence 
of a quasar-dominated UV/X-ray background (Haardt & 
Madau 1996; Madau & Efstathiou 1999); in this case the 
medium is highly photoionized, collisions are not impor- 
tant in determining the ionization structure (although we 
include them self-consistently) , and collisional excitation 
cooling of H i and He n becomes ineffective (Efstathiou 
1992). The cooling timescale for optically thin intergalac- 
tic gas at different overdensities (S = PbtPb = 1,5, and 10) 
and redshift z = 3.5 is shown in Figure [p it is longer than 
or at best comparable to the expansion timescale, 



is shorter than the expansion timescale. Here m e is the 
electron mass, ot the Thomson cross-section, a the radia- 
tion constant, and Tcmb = 2.725 (1+z) K the temperature 
of the CMB (Mather et al. 1999). Pregalactic outflows will 
propagate with typical velocities of a few tens of kms -1 
into a dense IGM which has been pre-photoionized by the 
same massive stars which later explode as SNe, and the 
expansion of the metal-enriched bubbles will be halted by 
the external pressure. By z = 3 any residual peculiar ve- 
locity would have been redshifted away by a factor of 2-3, 
the Lya forest would be hydrodynamically 'cold', and the 
intergalactic baryons would have relaxed again under the 
influence of dark matter gravity. 

(3) In a CDM universe, structure formation is a hierarchi- 
cal process in which non linear, massive structures grow 
via the merger of smaller initial units. Large numbers 
of low-mass galaxy halos are expected to form at early 
times in these popular cosmogonies, perhaps leading to an 
era of widespread pre-enrichment and preheating. The 
Press-Schechter (1974, hereafter PS) theory for the evolv- 
ing mass function of dark matter halos predicts a power- 
law dependence, dN/d\nm oc m^™ " -3 ^ 6 , where n e ff is the 
effective slope of the CDM power spectrum, n e s ss —2.5 on 
subgalactic scales. As hot, metal-enriched gas from SN- 
driven winds escapes its host halo, shocks the IGM, and 
eventually forms a blast wave, it sweeps a region of inter- 
galactic space which increases with the 3/5 power of the 
energy E injected into the IGM (in the adiabatic Sedov- 
Taylor phase). The total fractional volume or porosity, Q, 
filled by these 'metal bubbles' per unit explosive energy 
density E dN / d In m is then 

U oc m oc (dN/d\nm) ' oc m 

ll/M (8) 

Within this simple scenario it is the star-forming objects 
with the smallest masses which will arguably be the most 
efficient pollutant of the IGM on large scales. Note, how- 
ever, that since the cooling time of collisionally ionized 
high density gas in small halos at high redshifts is much 
shorter than the then Hubble time, virtually all baryons 
are predicted to sink to the centers of these halos in the 
absence of any countervailing effect (White & Rees 1978). 
Efficient feedback is then necessary in hierarchical cluster- 
ing scenarios to avoid this 'cooling catastrophe', i.e. to 

7 At the low metallicities (less than 1% solar) typical of Lya forest clouds, the thermal behavior can be modeled to a good approximation 
by a gas with primordial abundances (e.g. Sutherland & Dopita 1993). 

8 Assume, for example, that the chemical enrichment of intergalactic gas was due to the numerous population of Lyman-break galaxies 
(LBGs) observed at z = 3. With a comoving space density above + 1 = 25.5 of 0.013 h 3 Mpc - 3 (Steidel et al. 1999), a 1% filling factor 
would be obtained if each LBG produced a metal— enriched bubble of proper radius equal to about 140 ?t — 1 kpc. To fill such a bubble in 
5 X 10 8 yr, the ejecta would have to travel at an average speed close to 600 kms" 1 (for h = 0.5), with characteristic postshock temperatures 
in excess of 2 million degrees. 
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at all temperatures. 

While it is possible that some metals were dispersed in 
intergalactic space at late times, as hot pressurized bubbles 
of shocked wind and SN ejecta escaped the grasp of mas- 
sive galaxy halos and expanded, cooling adiabatically, into 
the surrounding medium, such a delayed epoch of galac- 
tic super-winds would have severely perturbed the IGM 
(since the kinetic energy of the ejecta is absorbed by in- 
tergalactic gas) , raising it to a higher adiabat and produc- 
ing variations of the baryons relative to the dark matter: 
Lya forest clouds would not then be expected to closely 
reflect gravitationally induced density fluctuations in the 
dark matter distribution, and the success of hydrodynami- 
cal simulations in matching the overall observed properties 
of Lya absorption systems would have to be largely coinci- 
dental. 8 In contrast, the observed narrow Doppler widths 
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prevent too many baryons from turning into stars as soon 
as the first levels of the hierarchy collapse. The required 
reduction of the stellar birthrate in halos with low circular 
velocities may result from the heating and expulsion of ma- 
terial due to OB stellar winds and repeated SN explosions 
from a burst of star formation. 

3. BASIC THEORY 

3.1. Dark matter halos 

To model the structural properties of subgalactic sys- 
tems we will neglect the gravitational potential due to vis- 
ible mass, and assume that virialized dark matter halos, 
formed through hierarchical clustering, have a universal 
(spherically averaged) NFW density profile, 



p(r) 



Pcrit 0~c 



cx(l + cx) 2 



(9) 



where x = r/r v i r , r vlr is the 'virial' radius of the sys- 
tem, i.e. the radius of the sphere encompassing a mean 
overdensity of 200, c is the halo concentration parameter, 
S c = (200/3)c 3 /F(c) is a characteristic overdensity, and 



F(t) = ln(l + t) - 



t 



1 + t 



(10) 



The mass of the halo within the virial radius is M = 
(47r/3)200p crit rA r . Equation (g) implies a circular velocity, 



2 GM(r)_ F(cx) 
c[ ' r c xF{cY 



(11) 



where V 2 = GM/r vir . Gas at radius r will escape from the 
gravitational potential well only if it has a velocity greater 
than 



v 2 e (r) = 2 



GM(r') 



dr' = IV, 



2 Hex) 



1+cx 



xF{c) 



(12) 



The escape speed is maximum at the center of the halo, 
v 2 (0) = 2V 2 c/F(c). 

To proceed further, we follow the algorithm described 
in the appendix of NFW and compute the concentration 
parameter (or, equivalently, the characteristic density con- 
trast S c ) of dark matter halos as a function of their mass 
in a SCDM model with fijif = 1, h = 0.5, as = 0.63, and 
n = 1. The algorithm assigns to each halo of mass M 
identified at redshift z a collapse redshift z co \\, defined as 
the time at which half of the mass of the halo was first 
contained in progenitors more massive than some fraction 
of the final mass. The assumption that the characteristic 
density of a halo is proportional to the critical density at 
the corresponding z co ii implies 



S C (M, z) cx 



1 + Zcoll 

1 + z 



(13) 



Since lower mass systems generally collapse at higher red- 
shift, when the mean density of the universe is higher, 
at any given time low-mass halos will be more centrally 
concentrated than high-mass ones. For M = 10 s h~ x M 
and z = 9, one finds (z co i\,c) — (12.2,4.8). At the same 
redshift, a 10 9 M© (10 7 M Q ) halo would have (z co n,c) = 
(11.9,4.7) [(z coll ,c) = (12.5,4.9)]. 



To study in detail the impact on the IGM of an episode 
of pregalactic star formation at 1 + z < 10, we will as- 
sume in the following a 'typical' concentration parameter 
of c = 4.8. At these epochs, the dark matter halo of a 
subgalactic system will be characterized by a virial radius 



0.76 kpc M 8 1/3 /T 



1 + z 
10 



a circular velocity at r v 



V. = 24 ms ; Ml 3 (±±1 

10 



1/2 



and a virial temperature 

_ GM »m p _ 4 5 2/3 J 1 + z 
vir " r vir 2fc U 8 ^ I 10 



(14) 



(15) 
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where fi is the mean molecular weight (p — 0.59 for a fully 
ionized hydrogen/helium gas) and Mg is the halo mass in 
units of 10 8 /i _1 M Q . The escape speed at the center is 



r. (01 = 77 kms" 1 U\ /Z 



1/2 



(17) 



Note that high-resolution N-body simulations by Bullock 
et al. (2001) indicates that high-redshift halos are actually 
less concentrated that expected from the NFW prediction. 
In this case we may be slightly overestimating the escape 
speed from subgalactic systems. 

3.2. Halo gas cooling 

If gas collapses and virializes along with the dark mat- 
ter perturbation to an isothermal distribution, it will be 
shock heated to the virial temperature and settle down to 
a density profile 



hiPgasM = hipo 



(imp 

2 JiTvir 



[v 2 e (Q)-v 2 e (r)} (18) 



(Makino, Sasaki, & Suto 1998). The central gas density 
Po is determined by the condition that the total baryonic 
mass fraction within the virial radius is equal to Jib ini- 
tially: 



Po 



200 „3 



c 6 Q b e J 



Pcrit f°(l + t) A /H 2 dt 



= 840 h~ 



(19) 



where A = 2c/F(c). At the virial radius /9 gas ( r vir) = 
0.00144po- 

Figures show the cooling function A and timescale 
icooi at the center and virial radius of an isothermal halo at 
z = 9, as a function of the halo virial temperature. Cooling 
by Lya line radiation becomes inefficient for gas tempera- 
tures below 2 x 10 4 K: in a gas of primordial composition 
at such low temperatures the main coolant is radiative de- 
excitation of the rotational and vibrational states of molec- 
ular hydrogen. Objects relying on H2 cooling are usually 
referred to as Population III objects. Primordial H2 is 
produced with a fractional abundance of /h 2 ~ 2 x 10~ 6 
at redshifts < 110 via the H~ formation channel. Starting 
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from this low value the H2 abundance increases in collaps- 
ing pregalactic clouds, molecular cooling becomes efficient, 
and stars can form. While most of the Lyman-continuum 
photons produced by these stars are quickly absorbed by 
the dense H 1 disk layers, radiation in the Lyman- Werner 
bands escapes into the IGM to form a soft UV cosmic back- 
ground which can photodissociate H 2 via the Solomon pro- 
cess, thus inhibiting further star formation. An estimate 
of the H2 equilibrium fraction under these conditions can 
be obtained by balancing the rates for the above processes: 

= JpL = 4. 8 x iQ-9 T v ° ir 88 (4^ Jlw.21)" 1 , (20) 

where Jlw,21 is the specific metagalactic flux 
in the 11.2-13.6 eV energy range in units of 
10~ 21 ergon -2 s -1 Hz -1 sr _1 , n e is the residual electron 
density after recombination, and we have taken the forma- 
tion (kf) and dissociation (kd) rates of Abel et al. (1997). 
Using the above relation, we have included the contri- 
bution to the cooling rate per particle due to molecular 
hydrogen (Martin, Schwarz, & Mundy 1996) in halos of 
different virial temperature. The derived cooling function 
and cooling timescale at z = 9 are depicted in Figures |||| 
for different values of Jim, 21- Note that here and below we 
are implicitly assuming that the universe is reionized at a 
redshift z < 9: this is because, after reionization occurs, 
there will be a universal background of photons above 
13.6 eV which inhibits the formation of dwarf galaxies 
both by reducing the cooling rate of gas within halos with 
?vir < 5 x 10 4 K, and also by suppressing the accretion of 
high entropy ionized gas into subgalactic fragments (e.g. 
Efstathiou 1992; Thoul & Weinberg 1996; Gnedin 2000). 

The fraction of the baryonic content of a halo that can 
actually cool and reach the center is determined by the bal- 
ance between the cooling and the dynamical timescales of 
the systems. As shown in Figure || for 4.3 < logXUr < 5.7 
rapid cooling by atomic hydrogen and ionized helium can 
occur at these epochs on timescales much shorter than the 
free-fall time, 

r dr' 
tff(r) = / = = 2.2 x 10 7 hT 1 yr 

Jo \/vl{r')-vl(r) 

x i}w) ^ J*dx>lF(cx')-F(cx)]-V 2 , (21) 

for a gas element at all radii r < r v j r , where T(cx) = 
[F(cx) + cx/(l + cx)]/[xF(c)]. Therefore, for masses in 
the range 10 s h~ x < M < 10 10 h~ x M Q , infalling gas never 
comes to hydrostatic equilibrium, but collapses to the cen- 
ter at the free-fall rate. Outside this mass range, i.e. when 
cooling is dominated by H2 (at the low-end) and free-free 
emission (at the high-end), the halo gas can be pressure- 
supported and form a quasi-static hot atmosphere. If we 
denote with r coo i the radius where the cooling time is equal 
to the free-fall time, a parameter /& can now be defined 
as the ratio between the gas mass within r coo i and the 
total baryonic mass within the virial radius, Q^M. This 
gas fraction is plotted in Figure |5| as a function of virial 
temperature. In halos with f b = 1 all the accreted gas 
can cool immediately, and the supply of cold gas for star 
formation is only limited by the infall rate. Conversely, 



in systems with /& <C 1, the supply of cold gas is regu- 
lated by the longer cooling timescale everywhere but for 
a small amount of gas in the very central region of the 
halo. When weighted with the steep PS mass function, it 
is the gas at the peak of the cooling curve - i.e. gas in sub- 
galactic systems with masses comparable to the masses of 
present-day dwarf ellipticals - than may be more readily 
available to be transformed into stars on short timescales, 
and give origin to explosive multi-SN events. 

In the adopted cosmology (SCDM with h = 0.5 and rms 
mass fluctuation normalized at present to erg = 0.63 on 
spheres of 8 hr 1 Mpc), M = 10 s h^ 1 M Q halos would be 
collapsing at z — 9 from 2-cr fluctuations. At this epoch, 
more massive halos with M = 10 10 h~ x M Q , while able to 
cool rapidly, would be collapsing from 3-er peaks and be 
too rare to produce significant amounts of heavy elements 
(in a gaussian theory the peaks above 3-<7 contain ~ 5% 
as much mass as those above 2-cr at a given epoch), un- 
less they were able somehow to form stars more efficiently 
than lower-mass objects. Even if this were the case, how- 
ever, the SN ejecta would escape the grasp of these more 
massive halos with less ease, and the fractional volume of 
the IGM filled by their metal bubbles be correspondingly 
small. Halos from 1-er fluctuations would be more numer- 
ous and contain most of the mass, but with virial temper- 
atures of only a few hundred degrees they would likely be 
unable to cool (fb <C 1) via H2 before reionization actually 
occurs. 

A detailed history of the chemical enrichment of the 
IGM should include the contribution to its metal content 
from all levels of the mass hierarchy at every epoch and 
is beyond the scope of this paper. Recent computations 
(Haiman et al. 2000; Ciardi et al. 2000; Ricotti, Gnedin 
& Shull 2001) trying to assess the importance of radiative 
feedback on Pop III halos have reached different conclu- 
sions. Although it is clear that a soft UV radiation field 
tends to prevent the cooling and collapse of these objects, 
and that star formation inside Pop III systems shuts off 
well before reionization, quantifying this effect is difficult 
as the field intensity depends on poorly unknown parame- 
ters as the star formation efficiency and the escape fraction 
of ionizing photons from the star formation sites into the 
IGM. Additional complications are introduced by locally 
produced Lyman- Werner photons which photodissociate 
H2 molecules in the parent halo gas (Omukai & Nishi 1999) 
and by a hard radiation components which could boost the 
formation rate of molecular hydrogen (Haiman et al. 2000) 
by increasing the electron fraction and thus feeding the H~ 
formation channel. Given these uncertainties we will take 
here a conservative approach and assume that all objects 
below the critical mass for Lya efficient cooling are sup- 
pressed. Below we will focus on the role played by what 
we have suggested might be the most efficient pollutant of 
the IGM on large scales, subgalactic systems with masses 
~ 10 s hr 1 M Q at redshift 9, when large numbers of them 
grow non-linear and collapse. 

3.3. Energy injection by SNe 

When gas cools well below its initial virial temperature 
and infalls to the center, it fragments into clouds and then 
into stars. Unless the stellar IMF is dramatically bottom- 
heavy, stars more massive than 8 M Q will form in the inner 
densest (self-gravitating) regions of the halo, eventually 
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releasing their binding energy in a supernova explosion, 
returning most of the metals to the ISM, and injecting 
about 10 51 erg per event in kinetic energy. We shall as- 
sume in the following that (i) all SN explosions take place 
at the center of the halo, and (ii) our stellar population 
forms in a time interval which is short compared to the life- 
time of massive stars. These two simplifying assumptions 
are delicate and need to be discussed. The first hypoth- 
esis implies spatial coherency among explosions, which in 
turn ensures that all the energy released by different SNe 
is used to drive the same 'superbubble', producing a cu- 
mulative, multi-SN event. This condition is not generally 
met in large galaxies like our own, where superbubbles oc- 
cur essentially at random in the disk. For objects with 
small circular velocities, however, the spatial coherence is 
essentially guaranteed by the fact that the size of the re- 
gion where star formation occurs is comparable or smaller 
than the characteristic scale of the bubbles. Indeed, if a 
rotationally-supported exponential disk with scale length 
rd forms in subgalactic fragments, and the specific angular 
momentum of the disk material is the same as that of the 
halo, then angular momentum conservation fixes the col- 
lapse factor to r v it/rd — V2/X, where A is the spin param- 
eter of the halo (sa 0.05 from N-body simulations, Barnes 
& Efstathiou 1987), and the equality applies to dark ha- 
los treated as singular isothermal spheres. Our fiducial 
M = 10 8 h' 1 M Q system at z = 9 has r vlr = 0.76 h^ 1 kpc 
and rd = 27 h^ 1 pc. If we assume, for simplicity, that the 
self-gravitating disk of mass Md follows an isothermal ver- 
tical profile with a thermal speed c s = 10 kms -1 , typical 
of gas which is continuosly photoheated by stars embed- 
ded within the disk itself, then its scale height at radius rd 
is h/r d = Vl6e X(M/M d )(c s /V c ) 2 , and the disk is rather 
thick (e.g. Wood & Loeb 2000). Within a scale height the 
gas density is approximately constant, and the radius R s 
of a supernova remnant in a uniform medium of density 
n and in the pressure-driven 'snowplough' stage is given 

by R s = {E^/n) 1 / 7 ^ 7 pc (McKee & Ostriker 1977). It 
will take then only 10 5 yr for an individual bubble to grow 
bigger than a disk scale length. 

The second assumption has to do with the ability of the 
collapsing system to transform the cold material into stars 
on timescales shorter than a few x 10 7 yr, the lifetime of a 
typical SN progenitor. A larger spread in the stellar birth 
times is likely to decrease the final number of SNe, as ion- 
izing photons from the first massive stars will reheat and 
ionize the infalling gas, thus inhibiting the formation of 
subsequent stars. It is difficult to assess the validity of such 
hypothesis, given our current poor understanding of star 
formation. While in a typical OB association in the Milky 
Way roughly three SN per million years will occur (Heiles 
1990), a much higher rate may be sustained in pregalactic 
systems due to the very short cooling timescale of the halo 
central regions. If we accept these two assumptions, which 
can only be met by subgalactic fragments, we can calcu- 
late the effects of energy deposition in a star-forming halo. 
Usually the distinction is made between blowout, a par- 
tial removal of the gas from a galaxy, and blow-away, in 
which the entire gas content is ejected back into the IGM 
(e.g. MacLow & Ferrara 1999). For a spherical system, 
as our admittedly idealized collapsing halo, the two terms 
are synonymous. Therefore, the superbubble will escape 



into the IGM only if it can lift out of the halo its entire 
gas mass. The amount of material transformed into stars 
can be parameterized as 

M* = n b f b f*M. (22) 

The cooled fraction of baryons, /{,, has already been dis- 
cussed above and found equal to unity in such halo mass 
range. There are no firm estimates for the star formation 
efficiency however, and we consider it as a free param- 
eter of our model. We will analyze three limiting scenar- 
ios: a low (/* = 1%), a medium (/* = 10%), and a high 
(/* = 50%) efficiency case. For h 2 n b = 0.019 and h = 0.5, 
M+ ranges between 1.5 and 75x 10 5 M . With a comoving 
space density of dark halos with masses above 10 8 h~ x M Q 
of about 80 h 3 Mpc~ 3 at z = 9 (see § 5) this corresponds to 
a stellar density parameter of 0.002 /*, i.e. between 0.4% 
and 20% of the total stellar mass inferred at the present 
epoch (Fukugita, Hogan, & Peebles 1998). We note that, 
on this assumption, /*, and hence the early luminosity of 
these systems, would exceed that which is predicted by 
the usual low-mass extrapolation of CDM galaxy forma- 
tion models (cf. White & Frenk 1991). In these scenar- 
ios /* is postulated to decline steeply in shallow potential 
wells, thereby reducing the population of low-luminosity 
galaxies and avoiding the so-called 'cooling catastrophe'. 
We also stress that the value /* = 50% is only meant to 
represent an extreme case and should really be considered 
as an upper limit to the star formation efficiency of our 
fiducial system. This is because, when most mechanical 
energy is injected by SNe after 3 x 10 7 yr (the main se- 
quence lifetime of a 8 M Q star) and SN-driven bubbles 
propagate into the galaxy halo quenching further star for- 
mation, only material within 0.4 r v i r has actually had time 
to cool, free-fall to the center, and form stars (Fig. [IJ). 
With the adopted density profile, the gas within the ra- 
dius where tg = 3 x 10 7 yr makes about 40% of the total 
baryonic mass QbM of the halo. In general, one would 
expect significantly lower efficiencies than 50% as the con- 
version of cold gas into stars is limited by the increasing 
fractional volume occupied by supernova remnants. Ad- 
ditional constraints on /* may come from the observed 
metallicity distribution of the most metal poor stars in 
the Milky Way halo, although the results will be subject 
to uncertainties in the low-mass end of the IMF. 

The IMF also determines the number of supernova 
events. Conservatively, we assume here a Salpeter IMF, 
with upper and lower mass cut-offs equal to M u = 120 M Q 
and Mi, respectively. The value of Mi is varied from 
the standard case Mi = 0.1 M Q to a value appropriate 
for a top-heavy IMF often predicted for very low metal- 
licity stars, Mi — 5 Mq, and up to the extreme case 
Mi = 30 M & , in which every star formed explodes as a 
supernova and many of them may eventually end their life 
as a black hole. 

3.4. Mechanical luminosity 

The halos under study produce a rather limited amount 
of very massive stars, and the IMF is then sampled in 
a stochastic manner. To determine the time-dependent 
mechanical luminosity injected by SN explosions when 
a fraction /* of the gas mass is converted into stars, 
we have repeatedly sampled the IMF by using a Monte 
Carlo method until the desired total stellar mass was 
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reached. This procedure yields for M; =0.1 M Q and /* = 
(1%, 10%, 50%) the following number of stars more mas- 
sive than 8 M Q (hence of SNe): TV = (1126, 11172, 55094), 
respectively. For this IMF, 1 SN occurs every v~ x w 
136 M Q of stars formed, with an average stellar mass 
{M s ) = 0.35 M©. For /* = 1%, the choice M z = 5 (30) M Q 
gives TV = 6076(3000), u' 1 = 24(52) M and (M s ) = 
13.5 (52.7) M Q . As we will see later, the differences in the 
efficiency and IMF will influence the fate and evolution of 
the ensuing superbubble and its metal content. 

To determine the main sequence lifetime, to^{M s ), of 
massive stars we have used a compilation of the data avail- 
able in the literature (Schallcr ct al. 1992; Vacca, Gar- 
many, & Shull 1996; Schaerer & de Koter 1997; Palla, 
private communication), and derived the approximate fit 



and energy conservation yield the relevant equations: 



Myr 



33 



\8 M @ J 



-3/2 



M s < 28.4 M , 
, -i/2 (23) 
3 - 4 Unfe) M S >28.4M . 



The extrapolation to masses larger than 60 M Q is quite 
uncertain; in general this massive stars are rare enough 
that this will not seriously affect our results. 

We can now derive the mechanical luminosity of the 
massive-star association driving the superbubble. Since 
all stars are assumed to be born coevally in a single burst 
of star formation, the spread in the SN energy deposi- 
tion is only due to the difference in toB for the var- 
ious masses. The mechanical luminosity is defined as 
L(i) = dE/dt, where E is the energy produced by the 
ensemble of SNe. We further assume that each super- 
nova releases Eq = 10 51 erg in kinetic energy. The deriva- 
tion naturally accounts for the stochastic behavior of L{t), 
which nevertheless has two clear features: i) a pronounced 
initial peak during the first 5 Myr after the burst, caused 
by the crowding of the explosions of the most massive stars 
which tend to have very similar ages (see expression for ton 
above); and ii) random oscillations around a mean value 
roughly equal to N E / uia,yi{toB{M s )} ■ 

3.5. Superbubble evolution 

In this and the following sections we will model the 
evolution of SN-driven bubbles as they blow out from 
our 10 s h' 1 M Q fiducial halo, allowing for radiative losses, 
gravity, external pressure, and thermal conduction. Cor- 
related multi-SN explosions will create large holes in the 
ISM of pregalactic systems, enlarging pre-existing ones 
due to winds from their progenitors stars. Most of the 
swept-up mass, both in the early adiabatic and in the fol- 
lowing radiative phases, is concentrated in a dense shell 
bounding the hot overpressurized interior, which yet con- 
tains enough mass to thermalize the energy input of the 
SNe. 

Superbubbles are canonically studied by using the thin 
shell approximation (Kompaneets 1960; Ostriker & Mc- 
Kee 1988), which has been checked against numerical sim- 
ulations giving excellent agreement (MacLow & McCray 
1988). The shell expansion, whose radius is denoted by 
R s , is driven by the internal energy, Eb, of the hot bub- 
ble gas. The pressure of such a gas (with adiabatic index 
7 = 5/3) is therefore Pb = Eb/2-KR^. Hence, momentum 



-{VspRs) = ^R 2 s {P b 



P) 



GM(Rs) 
R 2 



pV s . 



(24) 



^ = L{t) - ^R 2 s PbRs 



V s nUh{Tb), (25) 



where the subscripts s and b indicate shell and bubble 
quantitites, respectively. We have defined the volume en- 
closed by the shell as V s = (47r/3)i?f, the dots repre- 
sent time derivatives, and p is the pressure of the ambient 
medium, taken to be equal to the halo gas density within 
r v ; r , and to the IGM background density at z = 9 outside 
the virial radius. As at r v i r the halo is still about 60 times 
denser than the IGM, to avoid unphysical effects due to 
this jump we have allowed the two distributions to merge 
through an exponential transition of width A = 0.2r v i r . 
Finally, rig 6 A(T&) is the cooling rate per unit volume of 
the hot bubble gas, whose _average hydrogen density and 
temperature are tTh,& and Tb, respectively. The physical 
interpretation of the various terms is straightforward: in 
the momentum equation, the first term on the r.h.s. de- 
scribes the momentum gained by the shell from the SN 
shocked wind, while the second term corresponds to the 
momentum lost to the local gravitational field. The terms 
on the r.h.s. of the energy equation describe the mechani- 
cal energy input, the work done against the shell, and the 
energy losses due to radiation. 

It is important to remark here that we are actually 
neglecting the complicated web-like structure ubiquitous 
to three-dimensional cosmological hydrodynamical simu- 
lations (e.g. Cen et al. 1994; Zhang et al. 1995). In CDM 
universes, virialized systems form at the intersection of 
mildly overdense filaments, along which most of the mass 
accretion (inflow) actually occurs; outside the virial radius 
there will still be a power-law decrease in density within 
material that is flowing in, but of course the spherical as- 
sumption breaks down. Moreover, prior to the reioniza- 
tion epoch, a significant fraction of all baryons are not 
distributed in a low-density IGM but actually condense 
into numerous small halos with virial temperatures above 
the cosmological Jeans mass. Such 'minihalos' have not 
been yet resolved in large-scale cosmological simulations, 
and are expected to dominate the average gas clumping 
in the intergalactic medium (Haiman, Abel, & Madau 
2001). We have also ignored any effects due to possible 
inhomogeneities within the halo gas. The shallow slope of 
the CDM spectrum on mass scales < 10 8 M Q leads to all 
small-scale fluctuations going nonlinear almost simultane- 
ously in time. The evolution of these early halos will then 
be marked by frequent mergers, which could raise the gas 
to a higher adiabat by shock heating and yield complex 
velocity and density fields within the 'virial radius' (Abel 
et al. 2000). In this first assessment of pregalactic outflows, 
we shall assume for simplicity that much of their structure 
and hydrodynamics can be understood from spherical pro- 
files of the physical quantities. 

To determine the pressure P of the ambient medium we 
further assume that both the halo gas and the IGM are 
photoheated at a temperature of 10 4 K by the SN progen- 
itors. In general, the size of an intergalactic H n region 
around a galaxy halo will depend on the H-ionizing pho- 
ton luminosity N i} on the fraction / csc of these photons 
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that can actually escape the dense star formation regions 
into the IGM, on the IGM mean density njj, and on the 
volume-averaged recombination timescale, i roc . When the 
source lifetime t s is much less than (i rec , H^ 1 ), however, as 
expected for a subgalactic halo shining for a few x 10 7 yr 
before being blown away by SN explosions (or in the case 
of a short-lived QSO; Madau & Rees 2000), recombina- 
tions can be neglected and the evolution of the H n region 
can be decoupled from the Hubble expansion. The radius 
of the ionized zone is then 



Ri 



1/3 



47T71H 



(54 kpc) 



0.02 



-1/3 



10 



10 7 yr 



1/3 



(^ 52 /c S c) 1/3 , (26) 



where N, = 10 52 N 52 s" 1 is the ionizing photon luminos- 
ity due to 10 3 massive stars distributed according to a 
Salpeter IMF. Ri will then be larger than the final size of 
the SN-driven superbubble (derived below) for values of 
the escape fraction greater than a few percent. 

It is interesting to derive the timescale at which the 
postshock gas enters the radiative phase. One can show 
(Weaver et al. 1977) that when the ambient gas pressure, 
gravity, and cooling can be neglected, and L(t) = const, 
the solution of the above equations is 



R, = 



125 \ 1/5 
154tt / 



L^ 1 
P 



The temperature of the postshock gas is then 

2/5 



/ 3 iim p 



R 2 = 4 x 10 4 



L 



^Myr 71 



K. 



(27) 



(28) 



For /* = 1% and Mi = 0.1 M© one has L 38 = L/(10 : 
erg s _1 ) w 10; also, in the central region of our fiducial 
halo the total gas number density is n ss 10 cm~ 3 . At 
such temperatures and densities the cooling rate per par- 
ticle is rj 2 x 10~ 22 erg s" 1 . Hence the cooling time is 
icooi(£ = lMyr) « 1300 yr, much shorter than the dy- 
namical time of the system, ensuring that the shell forms 
very rapidly and justifying the use of the thin shell ap- 
proximation. These estimates also highlight the difference 
between a bubble produced by repeated SNe and a single 
point explosion. In the former case most of the energy of 
the bubble resides in the hot, very inefficiently radiating 
cavity gas generated by the continuous energy injection, 
while most of the mass is contained in a thin layer that 
collapses to form the dense, cool shell. 

The final step consists in the determination of the cool- 
ing rate of the bubble, which depends on the density and 
temperature of the hot interior. The thin-shell equations 
only provide a relation for the product of these two quan- 
tities (i.e. the pressure), so another physical relation must 
be derived. Most of the bubble gas mass comes from 
the conductive evaporation at the contact surface between 
the hot gas and the cold shell. The structure of conduc- 
tive/cooling fronts has been studied by several authors 
(e.g. Cowie & McKee 1977; McKee & Begelman 1990; Fer- 
rara & Shchckinov 1993). Under the conditions in which 



thermal conduction is unsaturated, the evaporative flow 
is steady, and radiative losses are negligible, the rate at 
which gas is injected from the shell into the cavity is 



dM c 
dt 



25fc b s 



C\T, 5/2 R S 



(29) 



where r\ — 6x 10 -7 (we have assumed a Coulomb logarithm 
equal to 30) is the classical Spitzer thermal conduction co- 
efficient. 

The interior distribution of temperature and density 
are known to obey a self-similar solution of the type 
f(x) = f c (l - C) q , where q = 2/5 (-2/5) if / repre- 
sents the temperature (density); f c is the value of the 
variable at the center and £ = r/R s , where r is the radial 
coordinate inside the bubble. According to this solution 
Pb oc nf,Tb = const, consistent with the hypothesis implic- 
itly made that SN shocks rapidly decay into sound waves 
inside R s . By integrating the density profile up to R s we 
obtain the total mass in the bubble, Mb = C^n^i? 3 , where 
C2 = (125/39)7r/^TOp. Differentiating this expression with 
respect to time and equating it to the evaporation rate 
(|29| ) , we finally obtain an equation for the evolution of the 
temperature Tb 



dTb 
dt 



3-^—R s 

Rs 



-P, 



23 Ci kT b m 
10 Ch R 2 S P S ' 



(30) 



This relation closes the system equations (p4[)-(|25|), and 
allows us to derive the temperature and density structure 
of the bubble and its cooling rate. 

4. NUMERICAL RESULTS 



The evolutionary equations (|2J), (|25|), Q30j) have been 
integrated numerically. We start by analyzing the low ef- 
ficiency case, /* = 1%, in some detail to highlight the 
relevant physics. As shown in Figure ^, the radius of 
the SN-driven bubble increases with time up to a final 
stalling value of Rf = 3.5 kpc, when pressure equilibrium 
is achieved. This happens in approximately 180 Myr, less 
than half of the then Hubble time; up to this point the evo- 
lution is then largely unaffected by the Hubble expansion. 
In the inital stages we find R s cx t 1A , a faster growth than 
given in (27), both because of the acceleration occurring in 
the halo density stratification and, less significantly, of the 
time-dependent mechanical luminosity. In the late stages 
of the evolution the shell evolves according to the 'snow- 
plough' (momentum conserving) solution, R s cx t 1 / 4 , as 
expected in the uniform IGM density field. Gas will be 
lost from the halo if its specific enthalpy exceeds its grav- 
itational binding energy per unit mass. The distinctive 
feature of blow-away can be seen in the velocity profile 
as a sudden jump of 10-15 km s _1 occurring at t = 40 
Myr. In all cases the velocity just before reacceleration 
is slightly lower than the escape velocity at the virial ra- 
dius, w e (r v i r ) = 48 km s _1 . The outer shock is radiative 
as long as it runs into the ISM of the protogalaxy. After 
breakout (t — 40 Myr) the inverse Compton cooling time 
is 230 Myr and the shock becomes adiabatic. (Compton 
cooling dominates since we assume the IGM to be pho- 
toionized, and collisional excitation cooling of H and He 
is accordingly suppressed.) For illustrative purposes, we 
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also show solutions where the gravity term was turned off, 
the pressure of the IGM was set to zero, and the halo 
gas mass was set to 50% of the standard value to mimic 
the possible presence of a cold galactic disk, which the 
(bipolar) outflow hardly couples to in numerical simula- 
tions (MacLow & Ferrara 1999). It is only in the last case 
that the final stalling radius differs significantly from our 
fiducial solution. Note that, except in the case of zero ex- 
ternal pressure, we actually stop the numerical integration 
of the equation at Mach number M. = 1 , as the shock de- 
cays into sound waves. We define the radius where that 
happens as the 'stalling radius', even if at that point the 
bubble pressure is still slightly higher than the external 
pressure of the IGM and the expansion will continue for a 
while longer. In practice, since the thin shell approxima- 
tion assumes a strong shock, our solution actually breaks 
down before A4 = 1. 

One feature of our model for the evolution of a wind- 
driven bubble is that conductivity erodes the inside of the 
cold shell and thereby mixes some swept-up material into 
the hot cavity. If the SN ejecta carry magnetic fields, then 
conductivity will be markedly reduced perpendicular to 
B and this process will be inhibited. To check whether 
our solution depends crucially on the value of the thermal 
conduction coefficient, we have run a case in which r\ in 
equation (30) was decreased by a factor of 100. We find 
that the temperature profile of the gas within the cavity is 
increased by less than 25%. Such a change leaves the rate 
of expansion of the cold shell essentially unaffected by the 
lower conductivity. This is because the (thermo)dynamics 
of the bubble is largely determined by the mechanical en- 
ergy input and adiabatic expansion. 

4.1. Varying the star formation efficiency 

We now explore the effect of varying the star formation 
efficiency /*. The results of the numerical calculations are 
shown in Figure ^. The evolution is qualitatively simi- 
lar in the three cases, but the final radius is significantly 
larger as /* increases: for example, the case /* = 10% pro- 
duces a bubble radius which is 3.3 times larger than for 
/* = 1%. Note that now the stalling point is reached only 
at ages comparable to the then Hubble time tjj, where our 
Newtonian solution breaks down. In the extreme case of 
/* = 50% the shell is still expanding at a velocity of 30 
km s" 1 at t = 0.5 in = 200 Myr. 

4.2. Varying the IMF 

Variations in the lower mass cut-off of the adopted 
Salpeter IMF will also affect our solution. Recall that 
the number of SNe produced by our three different choices 
of Mi = (0.1, 5, 30) M Q isftf = (1126, 6076, 3000), respec- 
tively, for /* = 1%; obviously, in the last case the pro- 
genitor stars are much more massive. The energy released 
into the ISM by these very high mass SNe may well be 
lower than assumed here, as part of the ejecta may fall 
back to form a black hole. Even if this is not the case, 
the ability of very massive SNe to drive outflows is lim- 
ited because of their short lifetime, as shown in Figure |[ 
The moderately top-heavy IMF case, Mi = 5 M Q , gen- 
erates a final radius of about 10 kpc, i.e. 2.5 times larger 
than for Mj =0.1 M Q . For the very top-heavy IMF with 
Mi = 30 Mq, the shell barely reaches 2 kpc, in spite of 
being driven by 2.7 times more SNe than in the fiducial 



scenario. The reason is that in this case the pulsed mass 
and energy input only lasts for ss 3 Myr, and the increased 
peak mechanical luminosity (L^g, = 10 3 ) is not enough to 
counterbalance the strong external pressure of the dense 
halo gas: the bubble enters the 'snowplough' phase well 
before blowing out into the IGM, and decelerates rapidly. 
To lift the halo gas out of the potential well it is then im- 
portant that the energy injection continues at least until 
blow-away occurs; after that, the shell will propagate to 
large distances while conserving momentum in the rarefied 
IGM. 

4.3. Bubble temperature 

The evolution of the gas temperature T c at the center 
of the hot, metal-enriched expanding bubble [recall that 
T(r) = T c (l — r/R s ) 2 / b inside the cavity] is shown in Fig- 
ure O. While a fraction of order 10% of the then Hubble 
time is spent at high temperatures (logT c = 6.6 — 7.2), 
after blow-away cooling sets in very rapidly. This behav- 
ior is common to all cases except the most extreme one 
with Mi = 30 M Q , which essentially lacks the first very 
hot phase. Bubbles that do not stall continue to expand 
and cool. Inverse Compton cooling is the dominant source 
of energy loss as can be inferred from the slope of the tem- 
perature profile, which is much steeper than expected from 
pure adiabatic expansion (in this phase R s oc i 1 / 4 , hence 
T c cx t- 1 / 2 if the gas is adiabatic). Thus, independently 
on whether a stalling radius is reached or not, all bubbles 
will continue to cool according essentially to the same law. 

5. DISCUSSION 

In the previous section we have modeled the evolution 
of SN-driven bubbles as they blow out from subgalactic 
halos and propagate into the intergalactic medium prior 
to the reionization epoch. We have seen that SN ejecta 
will escape the grasp of halos with virial temperatures 
2vir }t 10 4 " 3 K at z = 9 (when they collapse from 2-cr fluc- 
tuations) if a significant fraction of their baryonic mass 
can be converted into stars over a dynamical timescale. 
Depending on the star formation efficiency, IMF, and de- 
tailed physics of the expansion, we find that after about 
half of the then Hubble time these outflows will have pro- 
duced a warm, 10 5 K metal-enriched low-density region, 
with a size between a few and several ten times larger 
than the virial radius of the source halo, surrounded by a 
colder, dense shell which contains most of the swept-up 
mass. At t — 10 8 yr, typical shell bulk velocities range be- 
tween 12 kins" 1 (for /* = 1%) to 40 kms" 1 (/* = 10%) 
and up to 60 kms -1 (/+ = 50%). If the fraction of ioniz- 
ing photons which escape the dense sites of star formation 
into the intergalactic space is greater than a few percent, 
SN-driven bubbles will propagate into a pre-photoionized 
medium, and their expansion will be halted by external 
pressure, gravity, and radiative losses. When the veloc- 
ity of the outflow becomes subsonic, the enriched material 
gets dispersed by the random velocities in the IGM and 
the bubble looses its identity. 

Knowing the mass within the bubble and the metal yield 
per supernova we can now derive the metallicity volution 
of the bubble. We have used the compilation of Todini & 
Ferrara (2000) -essentially based on the Case A of Wooslcy 
& Weaver (1995) - who give the chemical composition of 



10 



IGM Metal Enrichment 



the ejecta for zero metallicity Type II SNe (at redshift 9 
the universe may be too young to host Type la SNe) as a 
function of their mass. At each time step of the simula- 
tion, we check for the massive stars that went supernovae; 
we sum over the mass of the heavy elements produced by 
them and divide by the mass of the hot gas within the 
bubble at that time. This is the typical metallicity inside 
the bubble if metals can be efficiently mixed within the 
cavity. Generally speaking, the metal-rich ejecta will be 
separated from the shell-evaporated gas by a contact dis- 
continuity. The main problem then is to disrupt this sur- 
face. When the cooling time in the cavity becomes shorter 
than its age, a cooled shell will form which rapidly becomes 
unstable (Gull 1973) to Rayleigh-Taylor (RT) and Kelvin- 
Hclmoltz (KH) instabilities. The 'fingers' produced by RT 
will penetrate deep into the high metallicity gas and may 
be eroded by KH instabilities due to the passage of rapidly 
moving hot gas, i.e. a mixing layer, perhaps creating a 
rather uniform distribution of metals inside the bubble. 

Further mixing, this time between the shell and the bub- 
ble gas, may occur after the expansion stalls. At that time, 
due to the restoring gravitational force due g on the shell 
from the DM halo, the acceleration vector points from the 
dense shell to the rarefied bubble, a RT unstable situation. 
The growth time of the RT instability on spatial scales A s 
is 



tRT 



^g (A - 1) 
A s (A + l) 



1-1/2 



-1/2 



(31) 



If the density contrast between the shell and the hot cav- 
ity gas is A = n s /rih ^S> 1, and g = 3 x 10~ n cm s~ 2 
at r = 10 kpc, £rt is less than 10 Myr for any reason- 
able thickness (tens of parsecs) of the shell. Hence mix- 
ing of shell material with the interior metal enriched gas 
may be quite rapid. Ferrara, Pettini, & Shchekinov (2000) 
have recently pointed out that additional processes related 
to mixing of the gas caused by time-varying gravitational 
tidal forces produced by halo interactions must be invoked 
if the presence of metals at lower redshifts is found to be 
ubiquitous. The time scale of the process is difficult to 
estimate as it depends on the peculariar velocities of the 
halos, which are predicted to be low at high redshift. Note 
that the metal enrichment of the hot gas inside the bubble 
will increase its cooling rate, which may become compara- 
ble with inverse Compton cooling. This is achieved around 
metallicity ss 0.5Z@ at z = 9. While our treatment may 
then underestimate the cooling rate, we have seen in the 
previous section that cooling has a weak effect on the over- 
all dynamics of the ejecta. 

The time evolution of the bubble average metallicity, 
Zb, is shown in Figure [l(] for the five cases studied. Its 
value is affected by the dynamical evolution of the bub- 
ble through the injection by evaporative of mass from the 
shell into the cavity. At later times, after blow-away, Zf, 
rapidly approaches a constant value determined by the 
fact that (i) metal production ceases with the last ex- 
plosion after t = 35 Myr, and (ii) the evaporative mass 
flow rate (i.e. the mass input in the bubble) drops, be- 
ing extremely sensitive (oc T 5 / 2 ) to the (rapidly decreas- 
ing) temperature. The metallicity becomes then the ratio 
of two constants, and is determined by the past dynam- 
ical evolution of the outflow. The final values are in the 



narrow range Z\, = 0.2 — 0.4 Zq, with the only excep- 
tion of the peculiar case with Mi — 30 M©, which has 
Zb — 3 x 10~ 2 Zq. In Table 1 we show the expected bub- 
ble and shell masses at stalling. The final ratio Mb/M s is 
of order 1%. Thus, if the heavy elements were completely 
mixed into the swept-up shell, then the resultant metal- 
licities would be lower than those derived above by about 
two orders of magnitude. It is possible, though, that the 
heavy elements may remain restricted to a region with a 
modest volume filling factor, and would not perhaps pen- 
etrate into the underdense medium between the halos. As 
already mentioned, the outcome depends on the onset of 
RT instabilities: non- linear development of these instabil- 
ities would lead to more widespread dispersal of the heavy 
elements through the IGM. The details are sensitive to the 
sound speed within the bubble when it stalls: this speed 
depends on how much small-scale internal mixing there 
has been within the bubble, due to conductivity, etc; but 
it could be substantially higher than the escape velocity. 
If so, then pressure-driven 'fingers' of hot enriched ma- 
terial (originating within the bubble) could propagate out 
into the IGM, readily reaching distances of order the inter- 
halo spacing if they maintained their identity and did not 
mix too rapidly. The effectiveness of mixing would depend 
on the extent to which conductivity is inhibited by mag- 
netic fields. It is worth noting that the bubble material 
(consisting substantially of supernova remnants) may well 
contain a magnetic field. If so, the dispersal of the first 
heavy elements and of the first ('seed') magnetic fields are 
intimately connected. If the 'fingers' are magnetised, then 
they may well propagate a long way before mixing with 
their surroundings. Whether the intergalactic heavy ele- 
ments are fully mixed is, however, an interesting question 
that is still open. The Lya forest observations imply that 
most lines of sight through each 'cloud' or filament en- 
counter some heavy elements. But it is not clear that the 
mixing is uniform. It is possible that, even at redshifts 
of 3-5, the IGM is inhomogeneous on very small scales. 
Even in the regions that have been enriched, the metals 
could be restricted to magnetised 'streaks' with a volume 
filling factor of as little as one percent. These would corre- 
sponding to the 'fingers' which we envisage having formed 
at z = 10, which would have been sheared and distorted 
by the subsequent gravitational clustering that led to the 
fully-developed forest at z = 3 — 5. 

What is the spatial extent of the ensemble of wind- 
driven ejecta? In the adopted flat cosmology {VIm = 1, 
h = 0.5, as = 0.63, n = 1, f2&/i 2 = 0.019) and according to 
the Press-Schechter formalism, the comoving abundance 
of collapsed dark halos with masses M ss 10 8 ft. _1 M at 
z = 9 is about 80/i 3 Mpc -3 , corresponding to a mean 
proper distance between neighboring halos of 15 /i -1 kpc 
and to a total mass density parameter of order 3 percent. 
In Figure O we plot the expected overall filling factor 
(or porosity) of metal-enriched material for different sce- 
narios, assuming a 'synchronized' (over timescales shorter 
than the then Hubble time) population of starbursting ha- 
los. As already mentioned, with a star formation efficiency 
of /* = 10% only a small fraction, about 4% percent, of the 
stars seen today would form at these early epochs. Still, 
the impact of such an era of pregalactic outflows could be 
quite significant, as the product of stellar nucleosynthesis 
would be distributed over distance that are comparable to 
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the mean proper distance between neighboring low-mass 
systems, i.e. filling factors > 20%. Larger filling factors 
may be obtained for larger efficiencies, moderately top- 
heavy IMFs, halos where a significant fraction of the gas 
is in a galactic disk and does not couple to the outflow (as 
matter is ejected perpendicularly to the disk), or from a 
population of more numerous sources - which would there- 
fore have to originate from lower amplitude peaks - but 
are not required by current observations. Ferrara, Pettini, 
& Shchekinov (2000) have recently considered the process 
of metal enrichment of the IGM, so it is useful to com- 
pare our main results with theirs. They determined the 
final radius of the bubbles by assuming that once blow- 
away has occurred, the shell enters the snowplough phase 
and it is finally confined by the IGM pressure. Essen- 
tially, their stalling radius is much smaller than derived 
in this work because they assumed a low cooling/star for- 
mation efficiency factor, faf* — 0.0025 (corresponding for 
a halo of total mass 2 x 10 s M to a star formation rate 
of 10 -3 MQyr -1 over 30 Myr), i.e. they studied outflows 
in a 'quiescent' star formation mode that is perhaps more 
suitable to larger systems, rather than to the starbursting 
subgalactic halos we have focused on in this work. 

It is interesting to note that, following an early era of 
SN-driven outflows, the enriched IGM may be heated on 
large scales to characteristic temperatures Tiqm ^ 10 4 6 K 
if the star formation efficiency is > 10%, so as to 'choke 
off' the collapse of further M < 3 X 10 8 h' 1 M Q systems by 
raising the Jeans mass. In this sense the process may be 
self-regulating. While more detailed calculations - which 
include the contribution to the IGM metal and heat con- 
tent from all levels of the mass hierarchy as a function 
of cosmic time - need to be done to fully assess the im- 
pact of subgalactic halos on the thermal (and ionization) 
state of the IGM, we argue here that it is unlikely that 
an early input of mechanical energy will be the dominant 
effect in determining the ionization state of the IGM on 
large scales (cf. Tegmark et al. 1993). We have shown in 
§3.5 that if the fraction of ionizing photons that escape the 
dense sites of star formation into the intergalactic space is 
greater than a few percent, pregalactic outflows will propa- 
gate into an IGM which has been pre-photoionized by the 
same massive stars which later explode as SNe. The rela- 
tive importance of photoionization versus shock ionization 



obviously depends on the efficiency with which radiation 
and mechanical energy actually escape into the intergalac- 
tic space. One can easily show, however, that, during the 
evolution of a 'typical' stellar population more energy is 
actually lost in ultraviolet radiation than in mechanical 
form (e.g. Madau 2000). This is because in nuclear burn- 
ing from zero to solar metallicity (Zq — 0.02), the energy 
radiated per baryon is 0.02 x 0.007 x mnc 2 ; about one 
third of it goes into H-ionizing photons. While the same 
massive stars that dominate the UV light also explode as 
SNe, the mass fraction radiated in photons above 1 Ryd 
is 4 x 10~ 5 , ten times higher than the fraction released in 
mechanical energy. Of course, a pre-photoionized universe 
at 10 4 K, heated up to a higher temperature by SN-driven 
winds, will recombine more slowly. The expected thermal 
history of expanding intergalactic gas at the mean density 
and with low metallicity is plotted in Figure [l^ as a func- 
tion of redshift for a number of illustrative cases. The code 
we use includes the relevant cooling and heating processes 
and follows the non-equilibrium evolution of hydrogen and 
helium ionic species in a cosmological context. The gas 
is allowed to interact with the CMB through Compton 
cooling and either with a time-dependent QSO ionizing 
background as computed by Haardt & Madau (1996) or 
with a time-independent metagalactic flux of intensity 
10~ 22 erg cm -2 s _1 Hz -1 sr~ 4 at 1 Ryd (and power-law 
spectrum with energy slope a = 1). The temperature of 
the medium at z = 9 - where we start our integration - has 
been either computed self-consistently from photoheating 
or fixed to be in the range 10 46 — 10 5 K expected from SN- 
driven bubbles with significant filling factors. The various 
curves show that the temperature of the IGM at z = 3 — 4 
will retain little memory of an early era of pregalactic out- 
flows. 
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Table 1 

Bubble and shell masses at stalling radius 





Mi/ Mq 


M b / Mq 


M S /M G 


M b /M s 


1% 


0.1 


1.0 x 10 5 


1.6 x 10 7 


0.7% 


10% 


0.1 


6.8 x 10 5 


5.5 x 10 7 


1.2% 


50% 


0.1 


2.5 x 10 6 


1.6 x 10 8 


1.5%+ 


1% 


5 


4.0 x 10 5 


3.2 x 10 7 


1.3% 


1% 


30 


1.0 x 10 5 


1.5 x 10 7 


0.7% 



Note. — f Evaluated at t = 200 Myr. 
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Fig. 1. — Cooling timescale t coo \ as a function of temperature Tigm for optically thin intergalactic gas at z = 3.5. The medium is assumed 
to have an overdensity of <5 = pb/~p b = 1 {dash-dotted line), 5 (solid line), and 10 (dashed line), primordial abundances, and to be irradiated 
by a quasar-dominated UV/X-ray background. 
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Fig. 2. — Equilibrium cooling rate at the center (solid lines) and virial radius (dashed line) of an isothermal halo at z = 9, as a function 
of virial temperature T v y, and in the absence of a photoionizing background (i.e. prior to the reionization epoch). The halo has an assumed 
baryonic mass fraction of fl^ = 0.019/i — 2 . The gas density dependence of the cooling function at high temperatures is due to Compton 
cooling off comic microwave background photons. The labeled curves extending to low temperatures show the contribution due to H2 for 
three assumed values of the metagalactic flux in the Lyman- Werner bands, 47tJlw ( m units of 10 -21 ergcm _2 s _1 Hz -1 ). 
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Fig. 3. — Equilibrium cooling time at the center and virial radius of an isothermal halo at z = 9, as a function of virial temperature T v 
and in the absence of a photoionizing background. Assumptions as in Fig. El with 47tJlw = 0.5 X 10 — 21 erg cm -2 s _1 Hz -1 . 




Fig. 4. — Halo cooling times at z = 9 as a function of radius for logT v j r = 4.3,4.7, 5.7 (dashed curves). The gas is assumed to be isothermal 
and in collisional ionization equilibrium. Dotted curves: same for halos with virial temperatures in the range 3.6 < logT v j r < 4.0, where 
cooling is dominated by H2 (47tJlw = 0.5 X 10 -21 erg cm -2 s —1 Hz -1 case). Solid curve: gravitational free— fall time of a NFW halo at z = 9 
(h = 0.5). 
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Fig. 5. — Fraction of the total halo gas mass at z = 9 that cools faster than the local free— fall time, i.e. never reaches hydrostatic 
equilibrium. In halos with 4.3 < logT v i r < 5.7 ft = 1, i.e. all the accreted gas can cool. For comparison, the fit for /(, derived by Machacek 
et al. (2000, MBA) from their hydrodynamics simulations is also shown. These numerical results have been obtained in the redshift range 
20 < z < 30, where because of the higher mean baryon density the rate of formation of H2 molecules is much faster than at the redshift of 9 
considered here. Some significant differences also arise from the different density and temperature profiles of Machacek et al. simulated halos. 
The gas in systems with T vn - < 5000 K (dotted vertical line) cannot cool faster than the local free-fall time at any radius. 
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Fig. 6. — Evolution of the shell radius and velocity in the fiducial case /* = 1%, Mi = 0.1 Mq (solid curves in both panels, repeated for 
clarity). All other curves show solutions with the same parameters when in turn: cooling is inhibited (small-dotted curves), the mechanical 
luminosity is kept at the constant value L38 = 10 (short-dashed curves), gravity is neglected (long-dashed curves), the IGM external pressure 
is set equal to zero (dot-dashed curves), and the halo gas mass is 50% of the standard value (large-dotted curves). A case with thermal 
conductivity suppressed by two orders of magnitude cannot be distinguished from the solid curve. 
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Fig. 7. — Same as Fig. ^j, but for the cases (f+,Mi) = (1%,0.1 Mq) (solid curves, repeated for comparison), (10%, 0.1 Mg) (dotted curves), 
and (50%, 0.1 Mq) (dashed curves). 
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Fig. 8. — Same as Fig. q, but for the cases (ft, Mi) = (1%,0.1 Mq) (solid curves, repeated for comparison), (1%,5 Mq) (dotted curves), 
and (1%,30 Mq) (dashed curves). 
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Fig. 9. — Evolution of the gas temperature at the center of the expanding bubble for the five cases (/*,Af;) = (1%,0.1 Mq) (solid curve), 
(1%, 5 Mq) (long-dashed curve), (1%,30 Mq) (dot-dashed curve), (10%, 0.1 Mq) (dotted curve), and (50%, 0.1 Mq) (short-dashed curve). 
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Fig. 10. — Average metallicity of the bubble Z^ (in units of solar) for the five cases shown in Fig. | as a function of time from the burst. 
If heavy elements get mixed up with shell material, the mean metallicity of the shell-bubble outflow will be two orders of magnitude lower 
than Z^. 
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Fig. 11. — Filling factor (porosity) of the metal enriched gas for the cases (1%,0.1 Mq) (solid curve), (10%, 0.1 Mq) (dotted curve), 
(10%, 0.1 Mq) where the halo gas mass is 50% of the standard value (large-dotted curves), and (50%, 0.1 Mq) (short-dashed curve). 




Fig. 12. — Thermal history of intergalactic gas at the mean density in an Einstein-de Sitter universe with Q b h 2 = 0.019 and h = 
0.5. Short dash-dotted line: temperature evolution when the only heating source is a constant ultraviolet (CUV) background of intensity 
10 — 22 erg cm - 2 s — 1 Hz - 1 sr _1 at 1 Ryd and power-law spectrum with energy slope a = 1. Long dash-dotted line: same for the time- 
dependent quasar ionizing background as computed by Haardt & Madau (1996; HM). Short dashed line: heating due to a CUV background 
but with an initial temperature of 4 X 10 4 K at z = 9 as expected from an early era of pregalactic outflows. Long dashed line: same but for a 
HM background. Solid line: heating due to a CUV background but with an initial temperature of 10 5 K at z = 9. Dotted line: same but for 
a HM background. 



